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<D Abstract 
c-j 

Under broad assumptions breaking of Lorentz invariance in gravitational theories 
^ leads to tension with unitarity because it allows for processes that apparently violate 

^ the second law of thermodynamics. The crucial ingredient of this argument is the 

ON existence of black hole solutions with the interior shielded from infinity by a causal 

^vq horizon. We study how the paradox can be resolved in the healthy extension of Hof ava 

<^ gravity. To this aim we analyze classical solutions describing large black holes in this 

theory with the emphasis on their causal structure. The notion of causality is subtle in 
this theory due to the presence of instantaneous interactions. Despite this fact, we find 
t ^ that within exact spherical symmetry a black hole solution contains a space-time region 

causally disconnected from infinity by a surface of finite area - the 'universal horizon'. 
We then consider small perturbations of arbitrary angular dependence in the black hole 
background. We argue that aspherical perturbations destabilize the universal horizon 
and, at non-linear level, turn it into a finite-area singularity. The causal structure 
of the region outside the singularity is trivial. If the higher-derivative terms in the 
equations of motion smear the singularity while preserving the trivial causal structure 
of the solutions, the thermodynamics paradox would be obviated. As a byproduct of 
our analysis we prove that the black holes do not have any non-standard long-range 
hair. We also comment on the relation with Einstein-aether theory, where the solutions 
with universal horizon appear to be stable. 
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1 Introduction and summary 



There are several reasons for the current interest in gravitational models with broken Lorentz 
invariance (LI). First, these theories may be relevant in the quest for consistent alternatives 
to general relativity (GR) that modify the laws of gravity at large distances. (This search is 
in turn prompted by the attempts to resolve the problems of dark matter and dark energy 
raised by cosmological data.) Second, the construction of these models provides a necessary 
framework for testing the nature of LI. These motivations are behind the Einstein-aether 
model [HE], where Lorentz breaking is realized by a unit time-like vector field, and the ghost 
condensation model [3], where LI is broken by a scalar field with time- dependent vacuum 
expectation value (VEV). More complicated patterns of Lorentz breaking are realized in 
models of massive gravity jl], see [5] for a review]^] 

More recently, this interest has been spurred by Horava's proposal [8J to construct a 
renormalizable model of quantum gravity by giving up LI. The main idea is that in the 
absence of LI the ultraviolet (UV) behavior of gravitational amplitudes can be substantially 
improved by the addition of terms with higher spatial derivatives to the action. This can be 
done while keeping the Lagrangian second order in time derivatives, thus avoiding problems 
with unstable degrees of freedom appearing in LI higher-derivative gravity [5] . The consistent 
implementation of this idea involves the notion of anisotropic (Lifshitz) scaling of the theory 
in the UV and indeed leads to a theory which is renormalizable by power-counting [8J. The 
renormalizability of the theory in the rigorous sense remains an open problem. 

Lorentz violating gravitational models generically contain new degrees of freedom besides 
the two polarizations of the graviton. These degrees of freedom survive down to the infrared 
where they can be conveniently described as a new Lorentz violating sector interacting with 
Einstein's general relativity. In this paper we will mainly focus on Hofava gravity. This 
model includes one new degree of freedom described by a scalar field <p(x) with non-zero 



Because of the latter property (p can be chosen as a time coordinate and thus acquires the 
physical meaning of universal time; we will refer to it as 'khronon'. The theory contains a 

*An important issue in theories with broken LI is to explain why Lorentz violation does not propagate 
into the Standard Model sector of particle physics where LI is tested with extreme accuracy. It has been 
proposed that LI may be protected by supersymmetry [S] . In [7] this protection mechanism has been realized 
in the context of the supersymmetric extension of the Einstein-aether model. 

2 We use the metric signature (+, — , — , — ). 



time-like gradient [T0| [TTj j^] 



dtf^O, {d^f > . 
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mass scale M* assumed to be just a few orders of magnitude below the Planck mass and 
which suppresses the operators with higher derivatives. These operators are not important 
for the description of the khronon dynamics and its interaction with gravity at low energy. 
The low-energy properties are captured by an effective 'khronometric' theory - a scalar- 
tensor theory satisfying certain symmetries [UJ. It is worth mentioning that, as it is always 
the case for effective theories, the khronometric model is in a sense more general than Hofava 
gravity: given the symmetries and field content, it can be derived independently as the theory 
with the smallest number of derivatives in the action. In this paper we concentrate on the 
'healthy' model proposed in [121 H3] which corresponds to the so called 'non-projectable' 
version of Hofava gravity. In this case the khronometric model possesses a symmetry under 
reparameterization of the khronon field, 



where / is an arbitrary monotonic function. Due to this symmetry the model differs essen- 
tially from ghost condensation, despite the fact that both are scalar-tensor theories with 
Lorentz violation)^} Instead, the model turns out to have a lot of similarities with the 
Einstein-aether theory [Til EL without, however, being completely equivalent to it. In- 
deed, in addition to the scalar mode the aether contains two more propagating degrees of 
freedom corresponding to the transverse polarization of a vector field. On the other hand, 
the khronometric model includes instantaneous interactions [11] that are absent in the case 



been analyzed in [TT| ITS] , where it was shown that for appropriate choice of parameters 
it satisfies the existing experimental constraint^] In [16] it was shown that extending the 
model by an additional scalar field with exact shift symmetry allows to naturally account 
for the cosmological dark energy. 

In this paper we study black hole (BH) solutions of the khronometric theory. These 
provide large distance (with respect to the scale M~ l ) solutions of the healthy Hofava gravity, 
and are expected to be produced by gravitational collapse. One expects BHs in Hofava 
gravity to differ from those of GR in several aspects. In GR a BH is characterized by the 

3 A model similar to the ghost condensation arises in the 'projectable" version of Hofava gravity [TUI lllj . 
However, in this case the khronon develops instability or strong coupling in the Minkowski background. 
4 This issue will be discussed in detail in Sec. [2] and Appendix |a| 

5 The emission of gravitational waves by binary systems containing sources with large self-energies has 
not yet been computed. However, for the systems where the radiation damping has been observed no big 
differences with respect to the computations of the weak fields regime are expected [TS] . 



V? ^ f(<p) , 



(2) 




The phenomenological consequences of the khronometric model have 
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existence of a horizon for light rays, which prevents the propagation of signals from its interior 
to the exterior. In this way, the horizon shields the central singularit}|^] of the BH from the 
exterior. However, in the presence of Lorentz violation the theory may contain excitations 
whose propagation velocity exceeds that of light and which thus can escape from inside GR 
horizons [T7] . Moreover, in Hof ava gravity the dispersion relations of the propagating degrees 
of freedom in the locally flat coordinate system have the form [8] 

where E and p are the energy of the particle and its spatial momentum, and c n , a n , b n are 
coefficients of order one depending on the particle species n. Stability at high momenta 
requires the coefficient b n to be positive. This implies that both phase and group velocities 
of particles indefinitely grow with energy and one might think that these modes can come 
from the immediate vicinity of the central singularity. Finally, as we discuss in Sec. |2j even 
within the low-energy description in terms of the khronometric model the theory contains 
a certain type of instantaneous interactions, which again might probe the BH interior down 
to the center. This suggests that Hofava gravity does not allow for BHs in the strict sense, 
characterized by the existence of regions causally disconnected from the asymptotic infinity. 
The purpose of this article is to clarify if this expectation is true. 

In our study we will be guided by the puzzles of BH thermodynamics arising in theories 
with Lorentz violation. As pointed out in [T8l IT9] for the examples of ghost condensation 
and Einstein-aether theories, one can construct gedanken experiments involving BHs that 
violate the second law of thermodynamics (see, however, [20] for a different point of view). 
Specifically, it is possible to set up processes that would decrease the entropy in the region 
outside the BH without any apparent change of the state of the BH itself. The second law of 
thermodynamics is intimately related to the unitarity of the underlying microscopic theory, 
see e.g. [2Tj . Thus its violation would constitute a serious problem, especially for a theory 
that, like Hofava gravity, aims at providing the microscopic description of quantum gravity. 
One can consider the following scenarios to recover the second law of thermodynamics: 

(i) The missing entropy is accumulated somewhere inside the BH. The first guess for the 
precise location of the entropy storage region is close to the central singularity. In 



6 It should be stressed that by singularity we always understand singularity from the point of view of the 
low-energy theory. In the full theory of quantum gravity the singularity is expected to be resolved by the 
effects that become important at high curvature (such as stringy effects in string theory or higher-derivative 
terms in Hofava gravity). 
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fact, this singularity is expected to be smeared off in the full theory at distances set 
by the microscopic scale and the entropy, in principle, can be stored in some high 
frequency modes localized in this smeared region. For this mechanism to work, an 
outside observer must be able to probe this region in order to make sure that the 
total entropy of the system does not decrease. This option seems plausible, given that 
Horava gravity contains arbitrarily fast modes including an instantaneous mode that 
persists at low energies. However, one must check that these modes can indeed escape 
form the center of the BH. In other words, one has to work out the causal structure of 
BHs in this theory. 

(ii) Another possible scenario to restore the thermodynamics assumes that BHs are not 
uniquely characterized by their mass, but instead can come out in many different 
configurations [21]. In other words, it supposes that BHs have a large number of static 
long-range hair. This hair would grow during the processes suggested in [181 EH] an d 
after measuring them an outer observer could decode the entropy fallen down into the 
BH. Notice that the difference with respect to the previous point is that the hair has 
a tail that can be measured outside the horizon. It has been demonstrated for several 
Lorentz violating theories including the Einstein-aether [22] , ghost condensate [23] and 
massive gravity |24j that spherically symmetric solutions with given mass are unique. 
As we are going to see, this is also the case for the khronometric model. This implies 
that the hair necessarily must be non-spherical. 

One of the purposes of this paper is to clarify which of these scenarios, if any, is realized 
in the healthy Hofava gravity. For this purpose, we first find spherical BH solutions in the 
khronometric model and analyze their causal structure. To simplify the analysis we will ne- 
glect the back-reaction of the khronon field on the metric. This approximation, valid when 
the dimensionless parameters of the khronon action are small, reduces the problem of finding 
the BH solution of the khronometric model to that of embedding the khronon field into a 
given metric backgrouncQ We show in passing that these solutions also describe spherical 
BHs in the Einstein-aether theory. Next we consider perturbations on top of these solutions. 
We demonstrate that no static long-range hair exists thus rejecting option (ii). On the other 
hand, option (i) is likely to work, though in a quite non-trivial manner. Within spherical 
symmetry we find that, despite the presence of arbitrarily fast, viz. instantaneous, interac- 
tions, the center of the BH is shielded by a causal horizon. This 'universal horizon' lies inside 

7 A similar approach was used in |24j for the analysis of BHs in massive gravity. 
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the Schwarzschild radius and its size linearly depends on the BH mass. All fields of the spher- 
ically symmetric solution are regular (analytic) at this horizon. However, in the spectrum 
of non- spherical time- dependent perturbations one finds certain modes with non-analytic 
structure at the universal horizon. These modes precisely correspond to the instantaneous 
interactions of the khronometric theory. Though at the linear order the universal horizon 
turns out to be stable under perturbations, we argue that the above non-analyticities will 
destabilize it at non-linear level, turning it into a finite-area singularity. One hopes that in 
the full Hof ava gravity this singularity is resolved into a high-curvature region of finite width 
accessible to the instantaneous and fast high-energy modes. In this way thermodynamics 
can be saved. 

It is worth stressing that the presence of instantaneous interaction is crucial for the type 
of instability discussed in this paper. Consequently, the universal horizon is expected to be 
stable in the Einstein-aether theory where all modes propagate with finite velocities. At 
present we are unable to suggest any resolution of the paradox with BH thermodynamics in 
this theory. 

Despite the vast literature on BHs in Hofava gravity, there has been only a few works 
dealing with the healthy extension in which we are interested. A class of spherically symmet- 
ric solutions of the healthy Hofava gravity was found in [25]. Those solutions differ, however, 
from the black holes considered in this paper. Ref. has obtained spherically symmetric 
BHs in the Einstein-aether and khronometric theories by numerically solving the coupled 
system of equations for the metric and aether (khronon) fielcj^J Our approach is different in 
that we simplify the setup to use analytic techniques whenever possible. This allows us to 
go beyond the spherical symmetry and study aspherical perturbations around BHs. Where 
our results overlap with those of they agree. 

The paper is organized as follows. In Sec. [2] we introduce the khronometric model and 
demonstrate that it contains instantaneous interactions. In Sec. [3] we find spherical BHs and 
show that those are also solutions for Einstein-aether theory. We also analyze their causal 
structure in this section. We turn to the analysis of the khronon perturbations about the BH 
in Sec. [4} Section [5] is devoted to discussion. In Appendix [A] we show that the instantaneous 
interaction is absent in the Einstein-aether model. Appendix [B] contains certain analytic 
results about spherical BHs. 

8 For previous studies of BHs in the Einstein-aether model see [STJ . 
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2 The khronometric model and instantaneous modes 



The khronometric action corresponding to the low-energy limit of Hofava gravity has the 
form [TTj . 



R + «(«%^) 2 + ^V/V,tt" + A(V/) 2 



(4) 



Here i? is the Ricci scalar and the unit vector is expressed in terms of the khronon field 

ip as, 

u^^=. (5) 

In Eq. Q M is a mass parameter related to the Planck mass and a, /3, A are dimensionless 
constant^} The previous action is the most general expression invariant under the transfor- 
mations ^ and containing only two derivatives of u^. Note that it formally coincides with 
the action of the Einstein-aether theoryp°| [TJ [2]. The difference, however, is that in the case 
of Einstein- aether the unit vector is treated as a fully dynamical variable, subject only 
to the constraint of been a unit time-like vector, and thus contains additional transverse 
degrees of freedom compared to the khronometric model where it is written in terms of the 
scalar field ip. This difference disappears for spherically symmetric solutions, since the aether 
vector is always hypersurface-orthogonal in that case [22] and thus can be expressed in the 
form ([5]). 

We assume that the parameters a, (3, A are small, 

a, f3, A < 1. (6) 



On one hand, this choice is motivated by the phenomenological bounds on the model (TTJ [15] . 
At the same time it drastically simplifies the rest of the analysis allowing to neglect the 
back-reaction of the khronon field on the metric. In fact, it is clear from Q that the energy- 
momentum tensor of the khronon field is proportional to the parameters a,/3,A . Thus its 
contribution into the Einstein equations is negligible when these parameters satisfy ([6]). This 
implies that the problem of finding a solution of the theory is reduced to solving the khronon 
equation of motion in an external (background) metric. 



9 The parameter A in (j4J) corresponds to A' in the notations of [11) . 

10 To be precise, the action of the Einstein-aether model contains an additional term proportional to 
V fl u l/ V fl u 1 ' . In the case of hypersurface-orthogonal aether ([5| this term is equal to a combination of the 
terms already present in the action pi). 
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A crucial property of the theory at hand is the presence of an instantaneous interaction 
[TT] . To understand its origin one can consider the propagation of small khronon perturba- 
tions in Minkowski space-time. It is straightforward to see that the Ansatz (p = t for the 
khronon background satisfies the equations of motion in this space-time. For the perturba- 
tions, we define a field x on top of this background, 

V = t + X , (7) 
for which one obtains the following quadratic action: 



a(S,x) 2 -(/3 + A)(A X -)'- 



Note that this action is fourth order in derivatives. However, the additional derivatives are 
purely space-like and the action describes a single wave mode with velocity 



As shown in [TT], this remains true also for curved backgrounds: the equation of motion for 
khronon perturbations always remain second order in time derivatives if the time coordinate 
is chosen to coincide with the background khronon field. 

Still the presence of extra spatial derivatives implies instantaneous propagation of signals. 
To see this, let us couple the khronon field to a source. The simplest source term that 
preserves the symmetries of the khronon field has the form 



s,„„« = j d v^sx , do) 

where S M is the source vector. To linear order around the flat background this reads 

SilLe = f d 4 x S% X , (11) 



which leads to the khronon exchange amplitude 

A = M- 2 S i (x)G ij (x - y)S*(y) , (12) 

where 



d p PiPj e 



-ipx 



G ' jU) 1 (27r)VM-(/? + A)p 2 ) ' (13) 
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where 9 is the step function and S the Dirac delta-function. A straightforward calculation 
yields the retarded propagator: 

gretfo) = d ( t ) 5 (\ x \- c x t ) x * x j , 0(t)0(\x\ ~ c x f ) ( 3x » x j ~ ^l x l 2 ^ t ( M ) 



4tt v /«(/3 + A) 



47TQ; 



Clearly, due to the second term the retarded Green's function does not vanish outside the 
khronon 'sound conej^]c x t = |x|. Moreover, the field extends to arbitrary spatial distance 
from the source immediately after the source is switched on, meaning that the signal prop- 
agates instantaneously. Note, however, that the instantaneous part of the signal builds up 
gradually starting from zero at t = and the maximum amplitude that is reached by the 
time t = |x|/c x decreases with the distance as l/|x| 2 . This is similar to the situation in 
Lorentz violating massive electrodynamics [28], [29] and massive gravity [211 EQ] • It is worth 
stressing that due to the existence of a preferred frame the instantaneous propagation of 
signals does not lead to any inconsistencies. 



Note that the instantaneous piece in ( 14 ) is traceless and transverse, 



Gf' ins \x) = <9 4 G[; Mns *(x) = . 



The latter property implies that the instantaneous contribution will vanish if the source 
S^(x) is a gradient of a localized scalar configuration. We also point out that, despite the 
similarity between the khronometric and Einstein-aether theories, the latter does not contain 



any instantaneous interactions. As shown in the Appendix [Aj the instantaneous piece in (14) 
is canceled by the contributions of the transverse modes. 



3 Spherical black holes 
3.1 Preliminaries 

According to the discussion of the previous section, our method to find BH solutions will 
consist in embedding the khronon field in an external BH metric. The latter satisfies the 
gravitational equations of motion in the vacuum. In the spherically symmetric case we will 
start with the standard Schwarzschild metric in the Schwarzschild coordinates, 

ds 2 = (l - -) dt 2 - - rW , (15) 

r 

11 Curiously, the Green's function does vanish inside the 'sound cone', as in the case of the ordinary massless 
field in four dimensions. 
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where 



dft 2 = d# 2 + sin 2 a ^ 2 



(16) 



is the area element of the 2-sphere and r s is the Schwarzschild radius. Using the Finkelstein 
coordinate 



v = t + r + r s In ( 1 



(17) 



the previous metric can be cast into the form regular at the Schwarzschild horizon, 

ds 2 = (l - -) dv 2 - 2dvdr - r 2 dtt 2 . (18) 

To simplify the khronon action, one observes that the curl of a hypersurface-orthogonal 
vector ([5]) identically vanishes, 



„ = e^ vXp u v V x u p = . 



This implies 



where 



Furthermore, 



2 



= = --u^u^ + {u v V v uJ 



(3 J d^v^V^W 4 = /3 J d 4 x^[(V^) 2 - R^uV] 



(19) 

(20) 
(21) 

(22) 



where R pv is the background Ricci tensor. Combining (20) and (22) we obtain that in the 



Schwarzschild geometry the khronon action can be written in the equivalent form 

+ A,„ ,^ 



a 



"(V/)' 



(23) 



where we have used that the Ricci tensor for the Schwarzschild metric vanishes. 



The khronon equation of motion, following from the variation of (23) with respect to ip, 
has the form of a current conservation, 



V„J" = , 



with 



J" 



X Su v 



(24) 
(25) 
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Here X = <9 M <^VV and 

J* = 8»- u»u v (26) 
is the projector on the hypersurface orthogonal to this implies 

= (27) 



We are interested in static, spherically symmetric, asymptotically flat solutions of (24). In 
the present context these notions require some clarification. The condition of staticity does 
not mean that the khronon field is independent of time (indeed, we saw that in the case of 
the Minkowski background tp grows linearly with time). Rather it implies that all quantities 
invariant under the symmetry ^ must be constant in time. In particular, this applies to 
the components of the vector u^. Combined with the requirement of spherical symmetry, 
this implies that the only non- vanishing components of the vector are ut, u r and that they 
depend only on the radial coordinate r. Note that we have chosen the position of the index 
T - down and V - up so that the corresponding components of the vector coincide in the 
Schwarzschild and Finkelstein frames: 

^ r I r I 

^ v i H | Finkelstein ^ I Schwarzschild 

For the khronon field itself, staticity and spherical symmetry imply that, up to a reparame- 
terization of the form ([2]), it can be cast into the form 

f = t + f(r) , (28) 

where f(r) is a function of the radius. Concerning the condition of the asymptotic flatness, 
we will impose, in addition to the flatness of the metric at r — > oo, the requirement that the 
khronon field tends to the same form as in Minkowski space-time, (p — t. This implies the 
boundary conditions at infinity: 

u t ->■ 1 , u r ->■ , VX ->• 1 at r ->■ oo . (29) 



We now prove that asymptotically flat static spherically symmetric solutions of (24) 
satisfy the stronger equation 

= . (30) 



Taking the first integral of Eq. (24) we obtain 



JT = ^ > (31) 
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where G\ is an integration constant. On the other hand, the explicit expression for J r in 



terms of reads, 



r 



au r u t d u t -\ — d r u t 
r 



(P + \)u 2 t d r [ d r u r + 



r 



where the components of the vector satisfy the unit norm constraint 



u; 



u 



r\2 



r 



(32) 



(33) 



After substituting (31) into (32) and expanding the r.h.s. of (32) at large r, the solution of 
the resulting linear differential equation reads 



Co 



u 



+ C 3 r 



(34) 



2(/3 + A) 

where C*2, C3 are integration constants. We see that the only solution compatible with the 



boundary conditions (29) is obtained for the choice C\ 
Finally, from (27) one finds that J* is also zerc^ 



C3 = 0. This implies J r 



0. 



We point out that Eq. (30) is the same as the equation of motion for the Einstein-aether 



model. Indeed, the latter is obtained by the variation of the action (23) with respect to 
the vector w M imposing that the variations of the vector components preserve the constraint 
u^u^ = 1. This leads to the appearance of the projector in the expression (25) for 
ja» Together with the fact that spherically symmetric solutions are always hypersurface- 
orthogonal [H] this proves that static spherically symmetric asymptotically flat solutions in 
the khronometric and Einstein- aether theories are equivalent. 



3.2 Solutions 

Our strategy to find the khronon configuration corresponding to a BH will be as follows. We 
solve the equation J r = for the vector w M , and afterwards reconstruct the corresponding 
khronon field. Introducing the notations 

U = ut, V = u r , £ = ^, (35) 

r 



12 Another derivation of the equivalence between (24 1 and (|30| in the spherically symmetric case, which is 



also applicable for non-static configuration, is presented in However, it is based on the assumption that 
all leaves of the constant khronon field that foliate the space-time are regular and simply connected. This a 
priori assumption appears too restrictive in the case of BHs. 
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equations (30) and (33) take the form, 



El 
u 2 - 



'■V 



+ 



v 2 



24 



(36a) 
(36b) 



where prime denotes differentiation with respect to £. Note that in the new coordinate the 
position of the Schwarzschild horizon is £ = 1, the region £ < 1 lies outside this horizon, 



while £ > 1 corresponds to the black hole interior. Expressing V from 13 

V=-y/U 2 -l + £ 
we obtain a single equation for the U component: 

U" + 



(36b) 



(37) 



c 2 x U 



f/ 2 (i-4)-i + £ 



{u/) 2 + (UU' + l/2f , 2(U 2 -l + 



U 2 -l+^ 



+ 



. 



(38) 



One observes that the denominator in the second term vanishes at the point £ c where 



[l-c 2 )U 2 (Q = l-Z c 



(39) 



This is nothing but the equation determining the position of the sound horizon for the 
propagating khronon mode. Indeed, we saw in Sec. [2] that in Minkowski space-time with the 
vector aligned along the time direction the wave component of the khronon propagates 
along the rays |x| = c x t. Denoting the tangent vectors to these rays by n M we find that they 
satisfy 

(40) 



M \2 



This equation has covariant form and straightforwardly generalizes to curved backgrounds 
with generic vector where it describes the sound cone for the khronon waves. The defining 
property of the sound horizon is that inside it all future-directed khronon rays point towards 



the center, n r < for any satisfying (40) and n v > 0, while outside of it the sign of 



n r can be positive. This implies that at the sound horizon the sound cone contains the 
vector = (1, 0, 0, 0). Substituting this into (40) and using the metric (18) we obtain (39). 



Clearly, for the subluminal khronon, c x < 1, the sound horizon lies outside the Schwarzschild 



13 The choice of the minus sign in front of the square root corresponds to the configuration where the 
foliation is 'infalling', which is what one expects for the black hole configuration. The plus sign would be 
relevant for the white hole case. 
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horizon, £ c < 1, while for the superluminal case c x > 1, the situation is opposite, £ c > 1; for 
c v = 1 both horizons coincide. 



We are interested in khronon embeddings with a regular future sound horizon 14 that are 
expected to appear as the result of a gravitational collapse. The regularity of the solution at 



the sound horizon implies that the term within the square brackets in (38) must vanish at 



£ c . This imposes a relation between the first derivative U' and the function U itself at this 



point, providing an additional boundary condition. Using (39) to simplify the formulas we 
obtain, 

1 



2(1 - 



1 + c, 



Sc 



(41) 



where we have chosen the branch of the square root that gives a regular result at c x = 1. 



We solved Eq. (38) numerically imposing the boundary conditions (39), (41) and 



U(£ = 0) = 1. 



(42) 



This was done using a shooting procedure where the shooting parameter is taken to be the 
value of the field at the sound horizon, U(£ c ). This determines the positions £ c of the sound 



horizon and the first derivative of the field there through (39), (41). Given U(£ c ) and U'(£ c 



we integrate Eq. (38) from £ c towards £ = to find U(0). Iterating this procedure we find 



the value U(£ c ) that allows to satisfy (42). Finally, we extend the solution inside the sound 



horizon by solving Eq. (38) from £ c to £ — > oo. In this way we obtain a unique solution for 
every value of c x . 

The resulting profiles U(£) are plotted in Fig. [ljfor several values of the khronon sound 
speed c x . A few characteristics of the solutions are listed in Table [Tj An important property, 
that is immediately clear from Fig. [TJ is the existence of roots of the function U(£) (i.e. 
points where Z7(£) crosses zero). In Appendix [B] the existence of at least one root is proven 
analytically for any value of c x . Moreover, the numerical analysis indicates that the number 
of roots is actually infinite with the function £/(£) exhibiting an oscillatory behavior around 
the zero axis. (The beginning of the oscillations is clearly visible in Fig. [T] for the curve 

14 Let us mention that in the superluminal case c x > 1 one can also find solutions that are regular 
everywhere except at r — (£ = oo) and do not possess the sound horizon at all. However, the numerical 
study of spherically symmetric collapse in the Einstein-aether model indicates that it always leads to the 
formation of a sound horizon [27] . Given the equivalence of the two theories in the spherically symmetric 
setting, one expects this to hold in the khronometric model as well. 

15 In the numerical implementation of the algorithm one has to set the initial data slightly away from £ c 

to avoid computational instabilities. 
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Figure 1: The time component U = Ut of the vector on the BH solution for several values 
of the khronon sound speed; the coordinate £ is related to the radius as £ = r s /r. 



corresponding to c x = 10; for smaller values of c x the oscillations lie outside the range shown 
in the figure.) This is in agreement with the findings of [26J. 

Let us denote the smallest root of U(£) by Its values for the numerical solutions at 
different c x are listed in Table. [Tj One observes that it mildly depends on the sound speed 
and satisfies > 1, > £ c . One can show that varies from 2 for c x — > to 4/3 when 
c x — > oo, see Appendix [Bj In the three-dimensional space the point corresponds to a 
two-sphere lying inside the Schwarzschild and sound horizons. We are going to see that this 
sphere plays a special role in the causal structure of the khronometric BHs. 

The last quantity presented in Table [l] is the derivative of the function U at = 
[/'(£*). These numerical data will be used below. 

3.3 The universal horizon 

Let us now study the khronon configuration corresponding to the BH solution. Taking the 
Ansatz 

<p = v + tp 1 (£) (43) 
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1 48 
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1 4zL 


n si q 


1 


1 


1.41 


-0 858 


1.5 


1.09 


1.38 


-0.906 


2 


1.14 


1.37 


-0.934 


10 


1.28 


1.337 


-1.032 


100 


1.328 


1.333 


-1.058 



Table 1: Several characteristics of the black hole solutions for different values of c x . See the 
text for definitions. 

and recalling the definitions (j5|, (35) we obtain the equation 

<*" = P(£-1)£W (44) 

where, to avoid proliferation of multiplicative factors r s , we have set r s = 1; this convention 
will be adopted from now on. It is straightforward to obtain the solution in the vicinity of 



£ = £*. Expanding the r.h.s. of (44) we obtain , 



Pi(fl« , 2r „ L= log (£,-£) at£«&- (45) 

We observe that the khronon field logarithmically diverges at ^. This implies that the 
leaves of constant khronon coming from the spatial infinity accumulate at £*. None of them 
penetrate inside the region £ > In other words, the leaves foliating the interior of the 
sphere £ = are disconnected from the spatial infinity. This fact is illustrated in Fig. [2] 
where the khronon foliation is superimposed on the part of the BH Penrose diagram covered 
by the Finkelstein coordinates. The situation is similar to what happens in GR for solutions 
containing Cauchy horizons. In both cases, the maximally extended solution contains more 
than one connected (asymptotic) region where boundary conditions should be specified. In 
the khronometric case, to know the solution at £ > £*, one needs to specify the boundary 
conditions for the instantaneous mode not only at io, but also at i + (see Fig. [2]). 

On the other hand, the sphere £ = simultaneously plays the role of the universal 
causal horizon. Indeed, the khronon field sets the global time in the model at hand. All 
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Figure 2: The leaves of constant khronon field (thin solid lines) superimposed on the upper 
half of the Penrose diagram of the Schwarzschild black hole. The thick solid line shows the 
universal horizon. 



signals, no matter how fast, can propagate only forward in this global time. In this way 
the configuration of the khronon determines the causal structure of space-time in Hofava 
gravity. From Fig. [2] it is clear that within this causal structure the inner region £ > lies 
in the future with respect to the outer part of the space-time. Thus no signal can escape 
from inside the surface £ = to infinity (null asymptotic region between i + and io) meaning 
that this surface is indeed a universal horizon, cf. [26J. 

It should be pointed out that within the spherically symmetric approximation that we 
have adopted so far the universal horizon is regular, despite the apparent singularity (45) 
of the khronon. Indeed, we have seen above that the field u^, which is the proper invariant 



observable of the theory, is smooth at £ = This implies that the singularity (45) can 
be removed by the symmetry transformation of the form (|2j). It is easy to see that the 
transformation 

^^^ = exp[(C 2 f/; V C 3 T) <p] 

does the job: the redefined khronon field is analytic at £*. However, in the next section 
we will argue that the universal horizon exhibits non-linear instability against aspherical 
perturbations of the khronon field, which turn it into a physical singularity. 
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4 Beyond spherical symmetry: perturbations 



4.1 Generalities 

The existence of the universal horizon seems to exclude the option (i) for the resolution of 
the thermodynamical paradoxes mentioned in the Introduction. However, this conclusion is 
premature: one has to analyze the stability of the causal structure depicted in Fig. [2] before 
making a definite statement. To this end we now study linear perturbations of the khronon 
on top of the BH solutions found in the last section. This will also allow us to explore the 
option (ii), namely possible existence of hair. As the background solution is not known in 



analytic form, we will characterize it by functions {/(£), V(£) as defined in (35) and obeying 



Eqs. (36) 



The analysis of the khronon perturbations must be performed in the preferred frame 
related to the background khronon foliation. In this frame the equations for the perturbations 
are second order in time, despite the presence of higher spatial derivatives [TOl [HI E] • Thus 
we introduce a new time coordinate r that coincides with the background khronon, 

T = v + <pi(£), (46) 



where (pi obeys Eq. (44). In the new coordinates the metric takes the form 
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Note that this metric is singular at £ = and thus only covers the region outside the 
universal horizon. This is sufficient for our purposes as we restrict to perturbations localized 
in this region. The complete khronon field is written as 

V = t + X (t,Z,6A) , (48) 

where x is a small perturbation. Due to the spherical symmetry of the problem, different 
spherical harmonics of x decouple from each other at the linear level, which allows to consider 
the equations separately for each multipole component xi- The latter obeys the relation 

A S 2 X i = -L 2 Xl , (49) 

where A52 is the Laplacian on the unit 2-sphere and L 2 = 1(1 + 1). To simplify notations we 
will omit the multipole label / on x in what follows. 



16 Recall that we work in the units with r„ = 1. 
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Instead of directly linearizing the equations of motion, we find it convenient to consider 



the quadratic action for the perturbations. Substituting (47), (48) into (23) after a straight 



forward (though somewhat lengthy) calculation one obtains 

L 2 



S {2) = 2nM 2 a J drd£ 



u 2 (x') 2 + ^x 2 - zeu'vx'x" - 2L 2 UV XX ' 
+ A(0(x'T + fliCOCxT + Q{i)x 2 



(50) 



The coefficient functions A, Bi, C\ are expressed in terms of U, V and their derivatives: 

A(0 =eu\V 2 - c\U 2 ) , (51a) 
= - 2£ A U 3 V 2 U" - £ 4 U 4 VV" - 2^U 2 V 2 U' 2 - 2£ 2 U 4 V 2 - 4£ 4 U 3 U'W (51b) 

- 4£ 3 U 4 W - 8£ 3 U 3 U'V 2 + L 2 i 2 U 2 V 2 - £ 4 U 5 U" - 2£ 4 U 3 V 2 U" 

- 4( - Q£ 4 U 4 U' 2 - 12£ 3 U 5 U' - 3£ 4 U 5 U" + 2L 2 i 2 U A + Q£ 2 U 4 V 2 - ^ 4 U 4 VV") , 
C,(£) = - L 2 [2£ UU'V 2 + i 2 UV 2 U" + 2£ 2 UVU'V + £ 2 <7 2 <y' 2 + t 2 U 3 U" (51c) 

- 4( - L 2 <7 2 + t 2 U 3 U" + 3^ 2 t/ 2 ^ /2 + MU 3 U' + 2U 4 - 2U 2 V 2 + £ 2 U 2 VV")] . 



Note that the coefficient A does not depend on the multipole number I. In deriving (50) 



and (51) we made use of the background equations of motion and integrated by parts to 



minimize the number of different x-structures appearing in the action. Note that the action 



(50) contains only two time-derivatives, as expected. 



Before analyzing the equation of motion following from (50) let us discuss the boundary 



conditions that must be imposed on the solutions. First, the field x must vanish at the spatial 
infinity, i.e. at £ — > 0. Second, x must be regular at the sound horizon £ = £ c . Note that 
the latter corresponds to a zero of the coefficient A(£). Determining the correct boundary 
conditions at the universal horizon is slightly more complicated due to the singularity of the 



metric (47) at To bypass this obstacle consider the linear perturbations of the vector u^. 

■ U) 



(1) 



eu 2 v x ' 



UK 



-eu 3 x ' 



u. 



(52) 



Here the positions of the upper and lower indices are chosen in such a way that the presented 



components are invariant under the coordinate change (46). Thus they also define the vector 
in the Finkelstein frame. We will require that the perturbations u^P in the latter frame 
are bounded at the universal horizon. This is compatible with the assumption that the black 
hole represents the end point of the gravitational collapse of a smooth initial configuration. 



From (52) we see that this requirement is equivalent to the condition that x diverges at 



19 



not faster than |£ — £*| 1 . We will see that imposing this condition is enough for our 
purposes. 



As a final remark we note that the expressions (51) considerably simplify in the case 



A =£ 4 (£ - 1)U 4 , (53a) 
Bi =2£ 4 (1 - £){U 2 + V 2 )V' 2 - 2£ 4 (1 - g)W' + 8£ 3 (1 - i)U 2 VV 



Q = -L 2 [{2ty 2 - 10££/ 2 + i 2 )UU' - 2t/ 4 + L 2 U 2 ] , 



(53b) 
(53c) 



where we have used the background equations (36) to simplify the final expressions 



4.2 Static perturbations: absence of hair 



Now we can prove that the khronometric BHs do not possess long range hah 17 At linear 
level such hair would manifest itself in the form of regular static perturbations of the khronon 
field (cf. [24]). These obey the equation 



(A(Ox"Y - (BtiWy + c,(Ox = o ■ 



(54) 



Let us count the number of free parameters of the general solution of Eq. ( 54 ) and compare 



it with the number of boundary conditions. Consider the asymptotics of the equation at 



spatial infinity, £ — > 0. Keeping the leading terms in the functions A,Bi,Ci, cf. (51), we 
obtain, 

- (£YT + 2L 2 (£Y)' + L 2 (2 - L 2 ) X = . (55) 

We will look for solutions of the power-law form, x ^ C" 5 - Substituting this Ansatz into (55) 
and solving the resulting algebraic equation for 5 one finds four roots: 



So 



-I , 6 3 = I + 1 , 5 A = -I - 2 . 



(56) 



Note that for I = we recover the asymptotics (34). For I > 1 the solutions corresponding 



to 82 4 grow at infinity and must be rejected. Thus we are left with a two-parameter family 



17 Of course, as in GR, BHs can have hair corresponding to angular momentum and possible gauge charges. 
We do not consider this standard hair, concentrating on those properties of BHs that are peculiar to the 
khronometric model. 
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of decaying solution; 



is 



X-Xi^ + Xs^ 1 ■ (57) 
These are further constrained by the boundary conditions at the universal and sound hori- 



zons. Expanding Eq. (54) at we obtain 

(£& - m\a - eovr + ((«* - i)gu: 2 (2£u: 2 - - ^ y)' = o . (58) 

Substituting the power-law Ansatz 

X oc (£ - £0 7 (59) 

we obtain for the exponent 

1 /l ^ / X 

71 = 0, 72 = -1, 7± = -o ±i/7 + 7^7772 ■ ( 60 ) 



Recalling that x must grow not faster than |£ — we obtain that the solution correspond- 
ing to 7_ is excluded. This gives one equation on the two parameters xu X3- One more 
equation follows from the requirement of regularity at the sound horizon £ = £ c . Thus in 
total we have two equations for two parameters. Assuming that this system is not degenerate 
we conclude that the unique solution is xi — X3 — implying the absence of hair. 

The fact that the equations following from the boundary conditions are indeed non- 
degenerate can be explicitly verified in the case of high multipoles, L 2 ^> 1. For the sake 
of the argument, we restrict to the case c x = 1, where the coefficients A, Bi, G\ simplify, see 



Eqs. (53). Keeping only the leading terms in L in Eq. (54) we obtain 



{eU\i - l)xT + L 2 ({-eU 2 V 2 + 2£ 2 f/ V)' - L 4 U 2 X = . (61) 

The form of this equation suggests to use the WKB method. Thus we search for the solutions 
using the expansion, 

x(0=exp[La(0 + ...]. (62) 



Substituting this into (61) and restricting to the leading order 0(L 4 ) we get, 



eU 2 (£ - 1)K) 4 + (-eV 2 + 2eU 2 ){a'f -1 = 0. (63) 

18 Thcrc is a subtlety in the case 1 = 1. The mode corresponding to Si is asymptotically constant and 
hence is sensitive to non-linear corrections. In principle, these corrections can make it diverge at £ — > 0. 
Whether this happens or not, is irrelevant for our argument. 
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This yields four solutions: 



\ l.± ^ «p ±L J — , (64) 

X2,± oc exp ± L J ^J^Tf: ■ ( 65 ) 

Note that xi,± are automatically regular at the sound horizon £ c = 1. We will see below 
that these solutions describe the instantaneous khronon mode in the black hole background. 
However, in the absence of sources they blow up either at spatial infinity or at the universal 
horizon (recall that U± < 0) in a way incompatible with the desired boundary conditions 
presented above. Thus they do not describe a valid static configuration. 

The two remaining solutions are at first sight singular at the sound horizon £ = 1. 
However, it is straightforward to check that the adiabaticity condition 



(7" 



<J 



A 2 



«1, 



(66) 



is violated at £ — >• 1, which means that the WKB approximation cannot be trusted in the 
vicinity of £ = 1. To work out the constraints imposed by the regularity at the sound 
horizon we have to solve Eq. (61) explicitly at £ ~ 1 and then match the solution to the 



WKB form (65) in a region where both approximations hold. In the vicinity of the sound 



horizon Eq. (61) takes the form, 

U A c {t - 1)X" + 2t/ c V + L'Utx" + L 2 (2Ut - Ul + AUlU' c ) X ' - L*U 2 cX = 



(67) 



where we have expanded all the coefficients at £ = 1 and introduced the notations, U c = U(l), 
U' c = U'(l). After introducing the new variable 

L 2 (£ 



C 



1 



Eq. (67) becomes to the leading order 0(L 6 ), 

A\ 2 d 3 X d^x 
^ dC 4 d£ 3 d£ 2 
The regular solution of this equation has the form 



X oc 



C 



(68) 



(69) 



(70) 



where J x is the Bessel function. The solutions (70) and (65) must be matched at L 2 |£| 
1. Using the asymptotics of J\ one finds that at £ < (£ < 1), the solution (70) contains 
both X2,+ and X2- components. The latter diverges at spatial infinity implying that this 
solution also cannot represent a static hair. 
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4.3 Time-dependent perturbations: non-analyticity at the univer- 
sal horizon 



We now return to the general case of time dependent perturbations. We will concentrate on 
large multipoles where it is possible to obtain approximate solutions using the WKB method. 
For clarity we will also restrict to the case c x — 1; the reader can easily get convinced that this 
restriction is not essential and that our qualitative results hold for arbitrary c x . Expanding 
the equation following from (50) to leading order in L and substituting x i n the form of a 
periodic function of time, 



we obtain, 

(£ 4 t/ 4 (£ - 1)XT + MH 2 U 3 V X "y + iu(eU 3 V X ')" - iuL\UVx)' - iuL 2 UV X ' 



(71) 



+ ([2L 2 eU 4 - L 2 eU 2 V 2 - cu 2 U 2 ] X ')' - (l a U 2 



u^L 2 



x = o 



(72) 



Upon substitution of the WKB Ansatz (62) this takes the form, 

2WVa' 
U 2 



i A U\i - l){a'f + 2int 2 U 3 V(a') 3 
-^U 2 V 2 -2eU 4 + Q 2 U 2 )(a') 2 

where we have introduced 1 ^ 



fl2 



1 



Remarkably, for any value of fl two of the solutions of Eq. ( 73 ) are 

1 



a 



±- 



(73) 



(74) 



(75) 



They coincide with the expressions obtained in the static case in the same WKB limit, 
cf. (64). This implies that the radial dependence of these modes is completely decoupled 
from their time-dependence, meaning that the field changes simultaneously all over the space. 
Thus these modes should be identified as mediating the instantaneous khronon interaction. 
If there are no external sources these modes are put to zero by the boundary conditions at 
spatial infinity and the universal horizon. We will discuss what happens in the presence of 
sources in a moment. 
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This corresponds to focusing on frequencies of O(L). 
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Before that let us consider the two remaining roots of Eq. (73) 



^ = a^ry U ± V 1 + — n* - J • (76) 

Recalling that V is negative and V(l) = —U(l) we see that the branch cr2,+ is regular at 
the sound horizon £ = 1, while o^,- has a logarithmic singularity. At infinity these modes 
correspond to the infalling and outgoing waves respectively, and thus the previous singularity 
just expresses the standard result that no outgoing modes can escape from the sound horizon. 
Both modes are regular at the universal horizon and cross it in the inward direction. Indeed, 
in the vicinity of we write, 



X2,± oc exp [ - iur + L cr 2) ±(£)] ~ exp 



log 

— 10J T 



e'™ , (77) 



where in the last equality we have switched to the Finkelstein coordinate v using the relations 



(46), (45). An implication of the above results is that the black hole is stable at the linear 
level with respect to the high-multipole perturbations. Indeed, an instability would manifest 
itself as a localized mode with positive imaginary part of the frequency. However, it is 
straightforward to check that for Imu > the solutions X2,± blow up either at £ = or 
£ = 1, implying the absence of localized modes. In this respect the situation is similar to 
the standard case of the Schwarzschild black hole in GR j2U [32], S3] . 

One might try to extend the stability analysis to arbitrary multipoles including the 
perturbations with I ~ 1. As we now discuss, this analysis appears unnecessary since there 
are strong indications that the universal horizon is anyway destabilized at the non-linear 
level due to the presence of the instantaneous interaction. To show this, let us look closer at 
the structure of the instantaneous mode. First, we notice that this mode can be separated 
from the part of the signal propagating with finite velocity for arbitrary multipoles provided 
that the frequency is large, u ^> 1. This is done by taking the limit u> — > oo in the 



equation following from (50) (alternatively, one can use Eq. (72)) and extracting the leading 
contribution. This gives the following equation for the radial dependence of the instantaneous 
mode: 

(U 2 X7 - ^ X = ■ (78) 
It is possible to show that this equation is equivalent to 
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where 

is the spatial Laplacian along the surfaces of constant khronon. Importantly, we again observe 
that the time dependence of the instantaneous signal completely factorizes from the spatial 



dependence. The general solution to Eq. (78) reproduces the asymptotics corresponding to 



5 2 , S3 at £ — > and 7+, 7_ at £ — > see Eqs. (56), (60). At high multipoles it reduces to 



the WKB solution given by ( 75 ) . 

As already noted, without any sources the instantaneous mode vanishes due to the bound- 
ary conditions. Let us study now what happens when there is a source for the khronon per- 
turbations outside the black hole. Physically, this can be either an external field interacting 
with the khronon, or the non-linear perturbations of the khronon left after gravitational col- 
lapse. To be concrete, we consider the situation when the source is located at a given distance 
from the center. Such a source produces instantaneous khronon perturbations that fall off 
at infinity and at the universal horizon. Due to the factorization property the perturbation 
in the vicinity of has the form, 

Xoc|£-£*l 7+ .f(<r), (79) 

where /(r) is the temporal profile of the source. The key feature of this expression is that 
it is non-analytic at £ = The singularity cannot be removed by a reparameterization of 
the khronon. The most straightforward way to see this is to consider the scalar invariants 
constructed from the vector u^. Let us take as an example 

/ = u^u*" , (80) 



where Ua V is defined in (21). At large frequencies, but still moderate multipoles, the leading 



linear contribution to this object is 

J (1) « A^UU'x' ■ (81) 

Though I^ 1 ' itself is finite (actually, vanishing) at £*, its derivatives of high enough order 
diverge. Indeed, consider the Zth derivative of along the trajectory tangential to the 
vector Its leading coordinate dependence in the vicinity of the universal horizon is 

( W A a A )'/Wocie-e,i 7+ "7 (0 (T). (82) 

It is easy to check using Eq. (60) and the data from the Table [I] that y+(t) < Z. In other 
words, we obtain that for a given multipole / the Zth and higher derivatives of the quantity 
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diverge. In particular, for the dipole I = 1 the divergence occurs already in the first 
derivative. It is natural to expect that these divergences will show up at non-linear orders of 
the perturbative expansion around the black hole background turning the universal horizon 
into a singular surface. 

Consider now a realistic dynamical collapse producing a BH. This system will never be 
perfectly spherical and will leave behind perturbations of the metric and the khronon present 
in the region outside the horizon and falling down as a certain power of time [32J. These 
perturbations will source the instantaneous mode which in turn will produce the divergences 
at the universal horizon. Thus we conclude that a realistic collapse, instead of producing a 
regular universal horizon, results in the formation of a finite area singularity at £ = £*. The 
Penrose diagram of the resulting configuration is depicted in Fig. [3] 




Figure 3: The Penrose diagram of the terminal configuration of a realistic collapse in the 
khronometric theory. The leaves of constant khronon field are shown by thin solid lines. 
Non-linear instabilities turn the universal horizon into a singular surface. The structure of 
the space-time inside this surface can be determined only within the full Horava gravity. 

The physical reason behind the instability of the universal horizon can be understood 
as follows. As already stated, for the instantaneous signal the universal horizon plays the 
role of a Cauchy horizon. For an observer falling down the BH and crossing the universal 
horizon the whole history of the universe outside the BH shrinks into a finite time interval. 
This leads to an infinite blue-shift of signals sent from the outside, these pile up at the 
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universal horizon and turn it into the singularity. The situation is similar to the instability 
of the inner Cauchy horizon in the case of Reissner-Nordstrom and Kerr black holes in GR 
[2H|35]. However, while in GR the signals destroying the Cauchy horizons are the standard 
waves propagating along light cones, in the khronometric case they must be instantaneous. 
Thus the presence of instantaneous interactions plays the crucial role in the destabilization 
of the universal horizon. Consequently, we expect this type of instability to be absent in the 
case of the Einstein-aether theory where all modes propagate with finite velocities. 

The singularity that replaces the universal horizon is naked in the sense that the com- 
plete determination of the khronon field requires imposing some boundary conditions on it. 
However, for those observers outside the black hole that are able to communicate only using 
finite-velocity signals (like those interacting only with the ordinary matter) the singularity is 
hidden by the corresponding sound horizons. Thus, despite the presence of the singularity, 
the evolution of the ordinary fields can still be unambiguously determined to the extent that 
their interaction with the khronon can be neglected. 

Finally, we do not know if the singularity discussed above is associated with a region 
of high space-time curvature, or if it is just an irregularity of the khronon foliation. One 
can conjecture that in the full Hofava gravity the singularity is resolved by the higher order 
terms, allowing for the complete determination of the dynamics. Admittedly, whether this 
happens or not is an open issue requiring further study. 

5 Discussion 

Let us summarize the picture that emerges form our analysis and speculate on its possible 
implications. It appears that in the healthy Hofava gravity a realistic BH formed as a 
result of a gravitational collapse has a central region - core - where the preferred foliation 
and possibly the metric are highly curved. The size of the core is proportional to the 
Schwarzschild radius with the proportionality coefficient mildly depending on the model 
parameters. The core appears singular from the point of view of the low-energy khronometric 
model. Optimistically, one may conjecture that the singularity is resolved once the full 
structure of the healthy Hofava gravity with the inclusion of higher-derivative terms is taken 
into account. 

The surface of the core can be probed from outside by the high energy modes and by 
the low-energy instantaneous interactions. In this sense the causal structure of the region 
outside the core is trivial: the causal horizons appearing for certain modes in the low-energy 
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approximation are absent once the whole spectrum of the theory is considered. Without 
invoking the full-fledged Hofava gravity it is impossible to tell whether the triviality of the 
causal structure will persist in the complete solution including the region inside the core. 
We believe that this is plausible, because the existence of the BH core accessible from the 
outside would eliminate the threat to BH thermodynamics raised by the violation of Lorentz 
invariance. Indeed, when considering the gedanken processes of jlHl [19] one would also have 
to take into account the changes in the entropy of the core: in the case of a trivial causal 
structure these changes are detectable from outside. Though at the moment we cannot 
estimate the entropy of the core one can speculate that its increase compensates for the 
decrease of the entropy in the part of the system outside the BH and in this way the second 
law of thermodynamics is saved. 

The existence of the BH core can also have profound consequences on the phenomenon 
of Hawking radiation. The (in) sensitivity of the latter to trans-Planckian physics is a long- 
standing issue, see e.g. [36j |37J |38l [391 HO]- Under broad conditions it has been shown 
that Hawking radiation is determined by the low-energy physics alone, in particular, by the 
properties of the sound horizons for low-frequency modes [SU [37J EH1 EH] • However, a key 
assumption in these works is that the quantum fields are in the vacuum with respect to the 
observers freely falling into the BH. This assumption is likely to be violated by the presence 
of the core just described: on the contrary, it seems more natural to assume that the fields 
are in the vacuum with respect to the rest frame of the core, which a priori is different from 
the free-falling one. It would be very interesting to work out what impact this can have on 
the spectrum of Hawking radiation. An extreme option would be that Hawking radiation 
gets completely suppressed. Hofava BHs would then behave as stable objects: a kind of 
'dark stars£3 

The drastic modification of Hawking radiation is also suggested by the following argu- 
ment. Consider adding to the healthy Hofava gravity a field carrying a global charge _ 
Imagine that a macroscopically large BH is formed and is endowed with large global charge. 
It is well-known that if the BH evaporates in the standard manner, the charge conservation 
will be grossly violated. The immediate way to see this is to note that the standard Hawk- 
ing radiation is charge-symmetric and thus particles and anti-particles are radiated in equal 
amounts leaving zero net charge after the BH evaporation. Even allowing for modifications 

20 This option does not exclude existence of a transient period just after the collapse when the BH would 
radiate nearly Hawking spectrum. 

21 We thank Sergei Dubovsky for suggesting to consider this setup. 
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of the radiation introducing a charge asymmetry will not save the day if the radiation re- 
mains (approximately) thermal with the temperature of the order of the standard Hawking 
temperature Tjj. Indeed, the BH starts giving away the charge only when Tjj exceeds m, 
the mass of the charged particles. Using the standard expression 

Tjj ~ M 2 p /M , 

one estimates the BH mass at this moment 

M ~ M 2 P /m . 

Conservation of energy then implies that the charge that can be given back by the BH is 
bounded above gH W2\ by 

Q < M 2 P /m 2 . 

On the other hand, the initial charge of the BH can be as large as 

Q ~ M /m , 

where M is the initial BH mass which, in its turn, can be arbitrarily large. Note that one 
can further relax the assumption of an approximate thermality of the radiation in the above 
argument. To show the non-conservation of charge it is sufficient to assume that (a) the BH 
evaporates and (b) that it cannot emit massive particles during most of its evolution. 

The charge non-conservation would prevent Hofava gravity from being a weakly coupled 
model of quantum gravity. Indeed, consider for concreteness the case when the charge is 
carried by a scalar field with a global U(l) symmetry. This symmetry is preserved at the 
perturbative level. Thus any violation of the charge conservation can stem only from non- 
perturbative gravitational effects. But in a weakly coupled theory, these are expected to 
be exponentially suppressed jl3] meaning that the charge must be conserved, at least with 
exponential accuracy. 

The only way to reconcile BH physics with the charge conservation, and thus save Hofava 
gravity as a candidate to a consistent quantum theory, is to admit that either even large BHs 
can emit massive charged particles, or that BHs do not evaporate at all. Clearly, both options 
would present qualitative departures from the standard picture of Hawking evaporation. 
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A Absence of instantaneous propagation in Einstein- 
aether theory 

We have seen in Sec. |2] that the khronometric model exhibits instantaneous interactions. 
From the mathematical perspective this is due to the fact that the equation describing the 
khronon is of fourth order because of two extra space-like derivatives. On the other hand, the 
Einstein-aether theory is described by second-order hyperbolic equations and one does not 
expect any instantaneous propagation in this case. This seems puzzling as khronon can be 
viewed as just the restriction of the aether to its longitudinal part. However, this is precisely 
the restriction that leads to instantaneous signals: in this Appendix we show that in the full 
aether theory the instantaneous piece is canceled by the transverse modes. 
The aether Lagrangian has the form 

M 2 r i 
C ae = -— [a(/V,«„) 2 + o^it,, W + ^V./V,/ + A(VX) 2 J • (83) 

This differs from the khronon Lagrangian Q by the presence of the term with the coefficient 
a'. The vector w M is subject to the unit-norm constraint 

«/ = 1 • (84) 

We are interested in the dynamics of perturbations in flat space-time around the background 
aether configuration 

uq = 1 , Ui — . (85) 



Due to the constraint (84) only space-like components of the aether are excited at the linear 



level. One separates the perturbations into the longitudinal and transverse parts, 

Ui — Ui = diX + uf , diuj- = . (86) 
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The longitudinal part is described by the quadratic khronon Lagrangian (|8j), where one 
should make the substitution: 



ol i — y ol -\- ol , 



The Lagrangian for the transverse component is 



a + a' . 



a fa ± 



(87a) 
(87b) 



(88) 



We now introduce coupling of the aether field to the source (10). The contribution of 



the longitudinal component to the exchange amplitude is given by Eqs. (12), (14) (again 



with the substitution (|87j)), and thus contains an instantaneous piece. However, a similar 
instantaneous contribution, but with the opposite sign, is present in the amplitude describing 
the exchange of the transverse mode: 
d 4 p Q-ipx 



ret , 



XI 



X) 



(2tt) 4 (a + a')pl - a'p 2 
9(t)5(\x\ - c±t) 



AiT\Ja'(a + cn')|x| 



PiPj 
P 2 



^(t)6'(|x| — c±t) f 3xjXj — Sij\x 
47r(a + a') 



?9) 



where c± 



is the velocity of the transverse mode. The last term in (89) cancels with 



that in (14) outside the wider sound cone, 



|x| = max (c_i_, c x ) t 



B Spherical solutions: analytic results 



In this appendix we present a few analytic results about the solutions of Eq. (38) (or equiv- 
alently, the system (36)) with the boundary condition (42) at £ = 0. 



We begin by proving that if there is a point £ c where (39) is satisfied, the solution will 



inevitably cross zero. In other words, all solutions having a sound horizon also possess a 



universal horizon. Let us assume the opposite. From (41) we observe that the derivative of 
U(£) at the sound horizon is negative, 

U'(Q < . (90) 
If U(£) stays positive at £ > £ c at least one of the following conditions must be satisfied: 
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(i) there is a point £1 > £ c where U' first changes sign, i.e. E/'(£i) = 0, U"(£i) > 0, 

(ii) £/(£) tends to a non-negative constant with U'(£) asymptotically approaching zero from 
below; this option implies that U" is positive at £ — > oo. 

We proceed to show that none of these two option can be realized. Consider the superluminal 
case c x > 1. Then it is easy to see that the combination 

F(£) = [/ 2 (£)(l-c x )-l + £ 

is positive at £ c < £ < £i- Indeed, this combination vanishes at the sound horizon, while its 
derivative 

F' = 2(1 - c 2 x )UU' + 1 

is positive in the above interval. According to the previous reasoning the factor in front of 



the square brackets in Eq. (38) is positive at the point £i, and it is straightforward to see 
that the quantity in the square brackets is positive as well (recall that E/'(£i) = 0). This 
implies that t/"(£i) is negative and we arrive at a contradiction. 



To exclude the option (ii) consider Eq. (38) at large £. For U" to be positive, the 



combination in the square brackets must be negative. This is possible only if (U') 2 > C 2 /£ , 
where C is some constant. This implies U' < — C/>/£ and thus U cannot asymptote to a 
non-negative constant. 

It remains to prove the statement in the sub luminal case c x < 1. The option (ii) is 
excluded by the same reasoning as above because the combination F(£) is positive at £ — > oo. 
However, the option (i) is trickier as it is not possible to argue that F(£) is positive in the 
interval £ c < £ < 1. Thus there may exist a turning point £i in this interval. To rule out 
this possibility we exploit the continuity of the solution in the parameter c x . If the universal 
horizon disappears for some values of c x , there must be a critical value c x such that £i exists 
for c x < and does not exist for c x > c x . For c x = c x the derivative £/'(£) touches zero 
at £i, 

U'fa) = C/"(£i) = . 



Then from Eq. (38) we obtain C/(£i) = 0. But this, combined with the requirement £i < 1, 



implies that the ^-component of the vector, Eq. (37), is undefined in the vicinity of this 



point. We have again arrived at a contradiction. This completes the proof. 
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The system (36) admits analytic solution in two limiting cases: 



and Cv — » oo. 



In the first case Eq. (36a) degenerates into U" = implying that the solution is a linear 
function, 

U = l-a£. 

The constant a is fixed by considering the first approximation in the small velocity c x < 1 



and imposing regularity at the sound horizon. In this limit (41) reduces to 



2£/'(£ c )£/(£ c 



-l. 



Solving for £ c from (39) in the same limit and substituting in the previous equation one gets 
a = 1/2. This yields the position of the universal horizon in this case, 



->• 2 at 



c x ->0. 



On the opposite extreme c x — > oo we can neglect the first term in (36a) and obtain the 
solution 

v = be, u 2 = \-i+b 2 e. (91) 

The solution exists at all £ only if b > bo, where bo = ^P. In the limit c x 



oo the condition 



(39) for the position of the sound horizon becomes 



U(Q = . 



For b > bo the functions £/(£) corresponding to (91 ) are strictly positive and thus the solutions 



do not possess a sound horizon. We believe that these solutions are unphysical and cannot 
be formed in the gravitational collapse. The physical solution that has a sound horizon 
corresponds tc 22 b = bo- The sound horizon in this case is located at £ c = 4/3 and coincides 



with the universal horizon. 
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